Experimental test of universal joint measurement uncertainty relations 
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The principle of complementarity is fundamental to quantum mechanics, and restricts the degrees 
of accuracy with which incompatible quantum observables can be jointly measured. Despite popular 
conception, the well-known Heisenberg uncertainty relation does not quantify this principle. Here 
we report the first experimental investigation of universally valid complementarity relations — also 
known as joint measurement uncertainty relations — and also derive a new, improved, complemen- 
tarity relation. Our tests exploited Einstein- Podolsky-Rosen-type correlations between two photonic 
qubits, to jointly measure incompatible observables for one of them. We demonstrate not only that 
the universally valid relations always hold, but that practical situations exist where the product of 
the measurement inaccuracies is sufficiently low to violate a widely used, but not universally valid, 
relation. 



Niels Bohr repeatedly emphasised that the fundamen- 
tal distinction between quantum and classical mechanics 
is the principle of complementarity, which states that the 
experimental arrangements required to accurately mea- 
sure two different observables are, in general, incompat- 
ible pp. Indeed, in replying to the famous critique by 
Einstein, Podolsky and Rosen on the completeness of 
the quantum theory [2], Bohr stated that "it is only 
the mutual exclusion of any two experimental procedures, 
permitting the unambiguous definition of complementary 
physical quantities, which provides room for new physical 
laws" [3]. 

It is commonly thought that Bohr's complementarity 
principle is captured by Heisenberg-type uncertainty re- 
lations such as AQAP > H/2, where AQ and AP are, 
respectively, the predicted spreads of perfect position and 
momentum measurements, arising from the underlying 
position and momentum spreads of a particle's wavefunc- 
tion. However, this relation does not reflect nor quantify 
the property that the mutual incompatibility of such per- 
fect measurements makes it impossible to carry them out 
at the same time. By contrast, a complementarity rela- 
tion should restrict the degree of accuracy with which 
the position and momentum can be measured simulta- 
neously. More generally, it should imply that no joint 
measurement procedure can generate perfectly accurate 
estimates of two incompatible observables. 

If two observables A and B are not jointly measur- 
able because they do not commute, one may still jointly 
measure estimates, ^4 e st ~ A and B cst ps B, where A est 
and -B cs t do commute. The key question then concerns 
the quality of these jointly-measurable estimates: how 
close are they to the original observables of interest? 
A natural measure of inaccuracy, for an estimate A est 
of some observable A, is the root mean square error 
e(A cst ) := ((A - A cst ) 2 ) 1/2 H 15] ■ This measure van- 
ishes for a perfect estimate, and more generally quanti- 
fies the degree to which the estimate statistically devi- 
ates from the quantity being estimated. It was shown by 
Arthurs and Kelly, in 1965, that e(Q ost ) e(P es t) > h/2 for 
any joint estimates of position and momentum which are 
globally unbiased, i.e., for which the average values of the 



estimates are equal to the average values of the corre- 
sponding observables for every state [BJ. Thus, the more 
accurately a globally unbiased measurement can estimate 
the position, the less accurately it can estimate the mo- 
mentum, and vice versa. The Arthurs-Kelly relation was 
subsequently generalised to the complementarity relation 

urn 

e(A cst )e{B cst )>^C, (1) 

for globally unbiased joint estimates of observables A and 
B. Here C := \([A, B})\, and is nonzero only when A and 
B are incompatible. Other complementarity relations — 
also called joint measurement uncertainty relations - 
have been given [51413] . but again are only valid for re- 
stricted classes of measurements. 

The abovementioned complementarity relations not 
only lack universal validity but, most importantly, they 
do not even apply in cases of considerable physical inter- 
est. For example, Einstein, Podolsky and Rosen (EPR) 
based their argument for the incompleteness of quantum 
theory on the case of two particles having very well de- 
fined relative positions and total momenta [2J. In this 
case, one can make very accurate joint estimates of the 
position and momentum of the first particle via, for ex- 
ample, a position measurement on the first particle and a 
momentum measurement on the second particle. Thus, 
the Arthurs-Kelly relation clearly fails in the EPR 
scenario. 

Indeed, it is only relatively recently that universally 
valid complementarity relations have been given, holding 
for arbitrary joint estimates of two observables in any 
measurement scenario. In particular, following related 
work by Ozawa on measurement-disturbance relations 
[14] . Hall showed that the inaccuracies of any estimates 
generated by a joint measurement of observables A and 
B must satisfy [T5] 

e(^ cst )e(B cst )+e(^ cst )AB cst +A^ cst e(£ ost ) > -C, (2) 

where AA 0St and AB cst denote the spreads of the esti- 
mates. Thus, the inaccuracies cannot both be zero for 
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incompatible observables A and B. This complementar- 
ity relation generalises ([!]) to arbitrary joint estimates, 
and is saturated in the EPR scenario [T5]. Hall also de- 
termined the form of the optimal estimates of A and B in 
any given measurement scenario, i.e., those functions of 
the measurement outcomes that yield the smallest pos- 
sible inaccuracies e(A est ) and e(B est ) [15] . Ozawa has 
given the formally similar relation |16j 

e(A est ) e(B cst ) + e{A cst ) AB + AA e(B cst ) >^C. (3) 

This differs from ([2]) in not depending on the spreads of 
the estimates per se, but on the spreads of A and B. The 
Hall relation is stronger than the Ozawa relation for the 
case of optimal estimates, since AA est < AA for such 
estimates [T5] , 

Here we provide the first experimental test of univer- 
sal complementarity relations, in an EPR-type scenario 
for which the Arthurs-Kelly relation is shown to be 
violated. Moreover, a new universal complementarity re- 
lation, which is stronger than both the Hall and Ozawa 
relations for optimal estimates, is proved here and simi- 
larly verified. 

A new universal complementarity relation. 

In addition to experimentally investigating comple- 
mentarity relations ([T])-([3]), we also tested the new re- 
lation 

Aj? cst + AB . AA cst + AA 1 
eldest ) 2 + e ( Scst ) 2 2 ' ( ' 

which we have derived and which is proved in ap- 
pendix A. Note that, unlike ([l])~([3]), it does not include 
an e(A cst )e(B cst ) term. Nevertheless, it similarly implies 
that e(A est ) and e(B est ) cannot both be zero for incom- 
patible A and B. 

Relation Q is typically stronger than both the Hall 
and Ozawa relations, and is always stronger in the case 
that optimal estimates of A and B are made, in any given 
measurement scenario, as shown in Appendix A. 

Measuring the inaccuracies. 

To test the validity of complementarity relations 
([4]) , it is necessary to experimentally determine the inac- 
curacies e(A est ) and e(-B es t)- This problem has been pre- 
viously discussed in context of measurement-disturbance 
uncertainty relations, which also involve the inaccuracy 
of an estimate. Ozawa proposed an approach in which 
e(A eBt ) is determined from the statistics of (separate) 
measurements of A and A est , on each of three different 
states of the system (defined by A and the state of in- 
terest) |17j . Lund and Wiseman proposed an approach 
not requiring preparation of additional states, in which 
a weak measurement of A is made prior to the measure- 
ment of A est |18j . allowing e(A cst ) to be determined from 
the corresponding weak- valued joint probability distribu- 
tion. These two proposals have been recently applied to 
the verification of measurement-disturbance uncertainty 
relations in Refs. [TS] and [2U] respectively. 



Our approach is a generalisation of the Lund- Wiseman 
proposal, in which the weak measurement interaction 
is replaced by one of arbitrary strength. To introduce 
the basic concept, consider two quantum observables 
K and L having respective eigenvalue decompositions 
K = J2k K kKk an d L = Yli^lLl- The corresponding 
'Margenau-Hill' joint quasiprobability distribution is de- 
fined by pmh(M) := {K k U + LiK k }/2 |Hj. It can take 
negative values for non-commuting observables, but is 
normalised and has marginal distributions correspond- 
ing to the probability distributions of K and L. The 
mean square deviation between K and L can be calcu- 
lated [IH1 [22] from the Margenau-Hill distribution via 



((K - L) 2 ) =J2( K k - Xi) 2 Pun(kJ). 



(5) 



k.i 



Further, if K is an observable prior to the application 
of an arbitrary quantum evolution (including measure- 
ment), and L an observable subsequent to its applica- 
tion, then the Margenau-Hill distribution itself can be 
determined experimentally via an initial weak measure- 
ment of the individual K k projectors postselected on a 
final measurement of L [TB]. Recalling that e(A est ) = 
{{A — Acst) 2 ) 1 / 2 , this allows the inaccuracies of joint esti- 
mates of A and B to be experimentally determined from 
weak measurements of the eigenprojectors of A and B 
(in separate experiments), postselected on the joint mea- 
surement. This approach also applies if L is replaced 
by any positive-operator- valued measure {L{\ having as- 
signed outcome values {A/} [15] , 

Our generalization of the Lund- Wiseman method re- 
laxes the requirement that the initial measurement inter- 
action be weak. Instead, it is only required that the nec- 
essary Margenau-Hill distributions can be obtained from 
the combined statistics of the initial and final measure- 
ments. This has the advantage of being simpler to imple- 
ment experimentally, although at the cost of increasing 
the inaccuracy of at least one of the joint estimates. Fur- 
ther, for our particular experiment, only one initial mea- 
surement rather than two weak measurements is needed 
to obtain e(A est ) and e(B cst ). The modification is in- 
spired by the contextual-value formalism of Dressel et 
at, which generalises many properties of weak measure- 
ments to arbitrary interaction strengths (2"31 [24] . 

In particular, suppose that in place of the initial weak 
measurement of the projector K k , the system undergoes 
an initial measurement X describable by a set of mea- 
surement operators {M m } [25], such that each K k can 
be written in the form K k = Y^ m a kmMmM m (for more 
general measurements, see Appendix B). The coefficients 
a km are called the contextual values of K k for the mea- 
surement context defined by X It follows that the 
joint probability distribution for outcome m from mea- 
surement X and / from the final measurement of L is 



p(m, I) = {MlLiM n 



(G) 
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and that the Margenau-Hill distribution for K and L is 
(k, l) = 5 L + hMlM m ). (7) 

m 

We have found that in many cases one can choose X such 
that the Margenau-Hill distribution ^ can be directly 
obtained either from the measured joint distribution ([6| 
or from (an assumed) theoretical description of L. For ex- 
ample, if [Li,M m ] = then p M n(k, I) = J2 m ctk m p{m, I). 
In all such cases the corresponding mean square deviation 
( (K — L) 2 ) can be evaluated via Eq. 
Experimental setup. 

To provide a discrete analogy of the EPR scenario dis- 
cussed in the introduction, we generated pairs of optical 
polarisation qubits having high fidelity with the entan- 
gled state 

\ip) = cos~/\HV) - sm~f\VH). (8) 

Here H and V refer to horizontal and vertical polarisa- 
tion, which are assigned as eigenstates of the Z operator 
|25j . For 7 = 7r/4 this state is maximally entangled, 
allowing perfectly accurate joint estimates of the X (di- 
agonal/antidiagonal) and Y (right/left-circular) polari- 
sations for the first qubit (via a measurement of Y on 
the first qubit and X on the perfectly correlated sec- 
ond qubit), analogous to the EPR example. However, 
this choice is not suitable for testing the complementar- 
ity relations since C = 2\(Z)\ = 2|cos27| and 
so trivially vanishes for 7 = 7r/4. Conversely, the choice 
7 = (or tt/2) maximises C, but corresponds to a fac- 
torisable state from which no information about the X 
polarisation of the first qubit can be gained via a mea- 
surement on the second qubit. Hence an intermediate 
choice, 7 = 7r/8, was made for our experiment. Quan- 
tum state tomography [23] on the state generated by a 
spontaneous parametric downconversion source gave a fi- 
delity of 97.4±0.3% with the desired entangled state (see 
Appendix C). 

The X and Y polarisations of the first qubit were 
jointly estimated using the outcomes of a direct measure- 
ment of Y on the first qubit and a measurement of some 
polarisation observable W of the second qubit (Fig. 1). 
To allow the corresponding inaccuracies to be determined 
by the method described above, the first qubit was also 
subjected to an initial semiweak measurement process X 
|27j . implemented using the polarisation dependence of 
Fresnel reflection from a glass slide [21] oriented at an 
angle fi= 75°. 

Horizontal and vertical polarisations were reflected 
with probabilities th and ry and transmitted with proba- 
bilities tff = 1—th andty = 1— ry (Fig. 1). Wave plates 
before and after the slide were adjusted to rotate the ba- 
sis of the semiweak measurement from Z to X, implying 
that the corresponding measurement operators M r and 
M t are Hermitian with M 2 = r H X + + r v X_ = 1 - Mf , 
where X± denotes |(1 ± X) and 1 is the unit operator. 
Values of r H = 0.1244±0.0001 and r v = 0.4645±0.0011 
were found by direct measurement. 




FIG. 1: Concept and implementation of the experiment. A 
source generates entangled pairs of optical qubits, and a joint 
measurement is used to simultaneously estimate the X and 
Y polarisations of the first qubit. The joint measurement 
comprises a semiweak measurement process X followed by 
a Y-polarisation measurement on the first qubit, and a W- 
polarisation measurement on the second qubit. The outcome 
of X is used to determine the inaccuracy, e(X cs t), of any es- 
timate X ost = f(W). The combination of the X and Y mea- 
surements denoted by y, is used to make a noisy estimate 
Y cst — y of the Y polarisation, with readily calculable in- 
accuracy e(Y eat ). Note this inaccuracy vanishes in the limit 
that A" is a perfectly weak measurement of X, with y = Y for 
this case. The process X uses a glass slide as a polarisation- 
dependent beamsplitter to make an semiweak measurement 
of X, and a flip mirror (FM) selects between the reflected 
and transmitted beams to enable the subsequent measure- 
ment of Y. Optical elements are coded by colour in the 
Figure, with representative elements labelled appropriately. 
QWP and HWP denote a quarter wave and half wave plate, 
respectively, FPC a fibre polarisation controller, and GT a 
Clan- Taylor prism. 

The contextual values — which represent the idea of 
eigenvalues generalised to arbitrary measurements [23] 
— are labelled a m , for x = ±1. These follow as a xm — 
I [1 ± x(mH + my)/ (rriH — my)} where + (— ) applies to 
the case m = t (to = r) (see previous section). Since the 
measurement process only acts on the first qubit, then 
[M m , W w ] = and the inaccuracy of any estimate of the 
form X cst = f(W) follows via Eqs. as 

epCst) 2 = ^2 [x - f(w)} 2 a xm p{m,y,w), (9) 

x,m,y,w 
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where p(m, y, w) is the joint probability of finding out- 
comes to, y and w for the binary measurements X , Y 
and W respectively. Hence, the inaccuracy can be di- 
rectly calculated from the measured reflectivities and 
joint probability distribution. 

The inaccuracy of the estimate Y est can be inferred by 
taking into account the effect of X on the first qubit under 
the assumption that the Y polarisation measurement de- 
picted in Fig. 1 is accurate. In particular, applying X be- 
fore the measurement of Y is equivalent to measuring the 
positive operator valued measure y = {Y y } on the initial 
state, with Y y = J2 m Mf n Y y M m [25]. This simplifies to 
Y y = k/2 + (1 - k)Y v , where k := 1 — ^r H r v - y/t H t v . 
Thus the effect of X is to add white noise of magni- 
tude k to the statistics of Y for the original state. Note 
that the value of n is a measure of the strength of X, 
with k = corresponding to the 'weak limit' rjj = ry, 
and k = 1 corresponding to a projective measurement 
of X. The measured values of rjj and ry above yield 
k = 0.0749 ± 0.008, moderately above the weak limit 
k — y 0. The Margenau-Hill distribution for Y and Y fol- 
lows as PMn(y, y) = 5 (YyYy + Y y Y y ) . The corresponding 
inaccuracy of the estimate Y cst = Y can immediately be 
calculated via Eq. (5) as 



e(Y ei 



^2[y - y\ 2 PMH{y,y) = 2k, 



(10) 



v-u 



irrespective of the input state. Note the inaccuracy 
vanishes in the weak limit k — > 0, as expected. 

Experimental results. 

The joint probability distribution p(rn, y, w) was mea- 
sured for several choices of the polarisation observable 
W, corresponding to the Bloch sphere angles 8 = 90° 
and (j) = 135°, 157.5°, 180°, 202.5°, and 225°. Note that 
9 = 90° and <f> = 180° corresponds to the observable of 
the second qubit that is most strongly correlated with 
the X polarisation of the first qubit, for state \ip) in ([8]), 
in the sense of maximising (X £g) W) . 

Two types of estimates for X, based on the outcome 
w = ±1 of the W measurement, were considered. The 
first type corresponded to simply estimating X to be w, 



X, 



simple 
est 



W. The second type of estimate exploited 



the tomographically-determined state of the source, p, to 
make the optimal estimate of the form X°£ t = f op t(W) 
with [T5] 



/optO) 



1 +wW 



1 +wW 



For example, for W = -X, i.e., 9 = 90° and (f> = 180°, 
this gave / opt (l) = 0.630 and / op t(-l) = -0.643. The 
corresponding inaccuracies were determined via Eq. {9), 
and are plotted in Fig. 2. It was also verified that 
the inaccuracy and spread of X°^ satisfied the Hall 
inaccuracy-dispersion relation |15j 
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FIG. 2: Inaccuracies of simple and optimal estimates o f X 
from a measurement of W, and verification of relation (111 
for the optimal estimate. Theoretical curves are given for the 
inaccuracies 

^simple) and e (x e °P t ) (upper and lower dotted 
lines), and for the spread AI C *' (dashed line). The solid line 
represents the tomographically determined value of AX, with 
the closely adjacent data points corresponding to the square 
root of the left hand side of Eq. (111. Error bars not shown 



are smaller than the size of the markers. 



for optimal estimates (Fig. 2). 

The inaccuracy of the estimate Y cst — Y was deter- 
mined via Eq. (To) as e(F est ) = (2k) 1 / 2 = 0.387 ± 002. 
Further, the spreads AX cst and AY^ st were calculated 
directly from the measurement outcomes, and AX = 
0.998 ± 0.002, AY = 0.9998 ± 0.0001, and C = 2\{Z)\ = 



1.422 ±0.008 from the tomographically-determined state. 

The above data allowed complementarity relations ([!])- 
Q to be tested, for the two types of estimate of X. 
The universal complementarity relations (J3) were val- 
idated and the Arthurs-Kelly relation violated (Fig. 3). 
The latter violation was expected for EPR scenarios, as 
discussed in the introduction. For the case of the opti- 
mal estimate of X (lower graph in Fig. 3), the results also 
verified the new complementarity relation Q is stronger 
than the Hall and Ozawa relations (|2) and ([3]), as ex- 
pected (see Appendix A). 

Discussion 

We have provided the first experimental verification 
of universal complementarity relations, that hold for any 
joint measurement of two observables (Fig. 3). This was 
done in an EPR-type scenario, for which the Arthurs- 
Kelly relation ([T]) is violated. We have, furthermore, de- 
rived and verified a new universal complementarity re- 
lation, Eq. Q, that is significantly stronger than the 
previously obtained relations ^ and ([3| for the case of 
optimal estimates, and have also verifed the inaccuracy- 
dispersion relation (111 for optimal estimates. 



(AX) 2 



(11) 



It would also be of interest to verify the universal com- 
plementarity relations for the case of the quadratures of 
a single-mode optical fields, for the continuously-valued 
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quantum dense coding [29 to the security of quantum 
cryptography |30j . and our work could have implications 
in all these areas. 
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FIG. 3: Experimental test of complementarity relations. The 
upper and lower graphs correspond to the simple and the op- 
timal estimates of X described in the text. Each graph shows, 
in descending order, the left hand sides of the Hall and Ozawa 
relations ([2| and Q (indistinguishable in the top graph); the 
new complementarity relation |4]);and the left hand side of the 
Arthurs-Kelly relation Q. The pink shading indicates the re- 
gion corresponding to violation of any of these relations. It is 
seen that the universal complemenarity relations |2|-([4| are 
satisfied whereas the Arthurs-Kelly relation Q is strongly vi- 
olated. Error bars not shown are smaller than the size of the 
markers. Solid curves are theoretical predictions. 



quadrature observables Q = a + a* and P(a — a')/i. 
This would allow the relations to be tested in the context 
more akin to the original scenario envisaged by Einstein, 
Podolsky and Rosen [5] , and have the advantages that (i) 
the righthand side of the relations is a fixed constant, in- 
dependent of the state, in this case; and (ii) the relations 
can be saturated for Gaussian states [15]. 

Our work represents an important advance in the 
quantitative understanding and experimental verification 
of complementarity, arguably the most important foun- 
dational principle of quantum mechanics. This principle 
underlies many aspects of quantum information technol- 
ogy, ranging from to entanglement verification [28] to 



Appendix A: Proof of Eq. Q 

Any two joint estimates of observables A and B may 
be represented by two commuting operators A ost and 
B cst on a suitable Hilbert space (via a Naimark exten- 
sion if necessary) [T51 [TBJ- It follows that 2[A,B] = 
[A - A cstl B + Best] + [A + A est ,B - B cst ], and hence, 
using the triangle inequality, that 

2\{[A,B}}\ = \{[A- A csU B] + [A- A cst ,B ost ] 
+[A,B-B est ] + [A est ,B-B est }}\ 
< \{[A-A est ,B])\ + \{[A-A est ,B est }}\ 
+ | ([A, B - Best]) | + \([A cst ,B - Bert]) | . 

Applying the Schwarz inequality (which lies behind all 
of these relations), \{[R,S})\ < 2 v /((B-r) 2 ) ((S-s) 2 ), 
with r and s suitably chosen from 0, (i?) and (S), leads 
to Eq. Q as desired. A corresponding measurement- 
disturbance uncertainty relation will be discussed else- 
where. 

Eq. Bl is always stronger than complementarity re- 
lations T2|) and (|3|), when optimal estimates A opt and 



B opt are made from the joint measurement results 
(corresponding to those functions of the measurement 
outcomes that yield the smallest possible inaccuracies 
e(A est ) and e(B cst )). This follows from the property 



(AA) 2 = (AA opt ) 2 +e(^ op t) 2 , (AB) 2 



for optimal estimates [T5] . 



(AB opt ) 2 +e(B opt ) 2 
(Al) 

In particular, AA opt < AA 



hence Eq. Q implies Eq. ^ 
her, defining a :— e(A opt )/AA 



and AB opt < AB, and 
for such estimates. Further, 
and 13 := e(B opt )/AB, then < a, ft < 1, and the 
difference of the left hand sides of relations (pi) and 
@ is e(A opt )ABh{ft) + AAe{B opt )h(a), where h[x) := 
(l/2)|Vl — x 2 — (1 — x)] is never less than zero. Hence 
Eq. Q also implies Eq. ([2| for optimal estimates, as 
claimed. 
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Appendix B: Inaccuracies and contextual values 

In our experiment, the initial measurement process X 
in Fig. 1 was describable by a set of measurement opera- 
tors {M m }, so that equations ^ and ffi were applicable 
to the calculation of the inaccuracies e(X C st) and e(Y" e st)- 
Here we indicate how this method can be extended to 
general X , via an extension of the definition of contex- 
tual values. 

In general, X will be described by some completely pos- 
itive trace-preserving map p — > x{p) = Em Xm{p)> where 
a measurement outcome m is obtained with probability 
Tr[xm(p)] and leaves the system in state Xm(p)/^[Xm(p)] 
|31jh Our experiment corresponded to the special case 
Xm(p) = MmpM^. Suppose that each projection oper- 
ator Kk for a Hermitian operator K can be written in 
the form K k = J^m^mXmW, where \ d m denotes the 
dual of Xm (i.e., Tv[A Xm (B)} = Tr[ X d m (A)B}), and 1 is 
the identity operator. The coefficients otum extend the 
definition of contextual values in |23j . and will again be 
referred to as the contextual values of Kk for the mea- 
surement context defined by X . Eqs. (|6| and ([7]), for the 
measurement joint probability and Margenau-Hill dis- 
tributions, then generalise to p(m,l) = (xfni^i)) and 
Pmh(M) = ^Y im a k m (x'L( 1 ) L i + LatW), respec- 
tively. As before, it is often the case that X may be 
chosen such that the Margenau-Hill distribution can be 
directly obtained either from the measured joint distri- 
bution ^ or from (an assumed) theoretical description 
of L. For example, if X operates on only one subsys- 
tem, and L is an observable of a second subsystem, then 
Xm and L may be rewritten in the forms Xm £5 I and 
1 <g> L, yielding p M n(k, I) = J2 m a kmP(m, I). In any such 
case the corresponding mean square deviation ( (K—L) 2 ) 



can be evaluated via Eq. ([5]), and hence the correspond- 
ing inaccuracies determined. 



Appendix C: Photon source and error analysis 

A 120 mW, linearly polarised continuous wave 404 nm 
laser diode was used to pump a pair of 0.5 mm thick 
Bismuth Borate (BiBO) crystals. The crystal pair is ori- 
ented such that the optical axes of the pair are crossed in 
what is known as a 'sandwich', and an additional BiBO 
crystal is inserted into the pump to precompensate the 
temporal walkoff between H and V polarisations in the 
downconversion crystal [321 133] . When pumped with di- 
agonally polarised light, spontaneous parametric down 
conversion leads to photon pairs at 808 nm in the maxi- 
mally entangled |$ + ) state. Counting in coincidence with 
a 3 ns window and 3 nm FWHM interference filters, pho- 
ton pair rates of ~ 2000 s _1 were observed. 

Uncertainties predominantly arise from poissonian 
counting statistics due to the random SPDC generation 
times, and slow, small thermal drifts in the fibre cou- 
pling leading to changes in rjj and ry during the data 
collection. The latter effect was characterized by measur- 
ing these quantites before and after the data collection, 
and observing the mean and spread of the values. Un- 
certainties in directly measured experimental values were 
combined using standard error propagation formulae to 
generate uncertainties in values calculated from the data. 
The small deviations between experiment and theory are 
attributed to a combination of inexact wave plate set- 
tings and other uncompensated small thermal drifts (e.g. 
fibre birefringence) in the experiment. 
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